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In the present study, optimization of rhombic drive mechanism used in a beta-type Stirling engine is 
performed based on a dimensionless theoretical model toward maximization of shaft work output. 
Displacements of the piston and the displacer with the rhombic drive mechanism and variations of 
volumes and pressure in the chambers of the engine are firstly expressed in dimensionless form. Sec¬ 
ondly, Schmidt analysis is incorporated with Senft’s shaft work theory to build a dimensionless ther¬ 
modynamic model, which is employed to yield the dimensionless shaft work. The dimensionless model 
is verified with experimental data. It is found that the relative error between the experimental and the 
theoretical data in dimensionless shaft work is lower than 5.2%. This model is also employed to inves¬ 
tigate the effects of the influential geometric parameters on the shaft work, and the optimization of these 
parameters is attempted. Eventually, design charts that help design the optimal geometry of the rhombic 
drive mechanism are presented in this report. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Stirling engines are potentially competitive when applied in a 
number of emerging engineering technologies since the engines 
are referred to as the external-combustion engines which are 
suitable for a variety of external heat sources. For example, the 
engines using nuclear energy could power deep-space exploration 
probes [ 1 ], using fossil or hydrogen fuel could drive a vehicle like an 
automobile [2], or a submarine [3], and using a solar dish collector 
could produce electricity [4]. A typical solar dish apparatus consists 
of a large dish aimed at the sun to reflect the rays into the focus 
point, which is located at a thermal receiver of the Stirling engine. 
The engine operates by converting the heat received from the solar 
dish into mechanical work. The work output of the Stirling cycle is 
then used to drive a generator and create electricity. Furthermore, 
the electrical power generated by the Stirling engine can be further 
used to electrolyze water for producing hydrogen. The hydrogen 
gas could be stored properly and then fed into a fuel cell for elec¬ 
trical power generation when needed [5], 

As reviewed by Senft [6], the Stirling engine was named after its 
inventor Robert Stirling who invented his first engine in 1816. The 
first successful theoretical analysis was presented by Schmidt [7] in 
1871. Schmidt’s analysis led to a simple closed-form solution for 
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indicated work of the engine. Lately, a large numbers of Stirling 
engine prototypes with various configurations, which have been 
categorized in to a-, p-, and y-type Stirling engines, have been 
created [8], In addition, some single-piston designs, such as 
thermal-lag Stirling engine [9] and pulse tube engine 10], are also 
presented. On the other hand, development in thermodynamic 
models [11—14], optimization methods 15—18], and measurement 
methods for engine performance [19—21] is also advanced by the 
progress in prototype building. Based on Schmidt’s model, a group 
of researchers, such as Finkelstein 15], Kirldey [16] and Walker [17] 
investigated the optimal phase angle and swept volume ratio for 
various types of Stirling machines. Senft [13] introduced a concept 
of effectiveness of mechanism and proposed a relation between the 
forced work and the shaft work of the engine. Later, the same 
author [18] combined the Schmidt theory and the shaft work the¬ 
ory to carry out the optimal combination of the phase angle and the 
swept volume ratio for the y-type Stirling engine. Recently, theo¬ 
retical analysis was extended to the irreversible cycle analysis. 
Cheng and Yu [12] proposed a thermodynamic model for the P-type 
Stirling engine with rhombic-drive mechanism by taking into ac¬ 
count the effectiveness of the regenerative channel as well as the 
thermal resistances of the heating and cooling heads. The ther¬ 
modynamic model was then incorporated with a dynamic model to 
perform a dynamic simulation model for the P-type Stirling engines 
with cam-drive mechanisms by the same group of authors [22], 
According to studies of and Kirldey [16] and Cheng and Yang 
[23 , it is recognized that among a-, P-, and y-type Stirling engines, 
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Nomenclature 

£ 

(R - l)lr 



6 

crank angle (deg) 

A 

cross section area (m 2 ) 

K 

swept volume ratio ( l/ sc /V se ) 

e 

I/r 

T 

temperature ratio (T c /T e ) 

E 

effectiveness of mechanism 

X 

total dead volume ratio 

l 

length of linkage (m) 



m 

mole number of working fluid (mol) 

Superscripts 

p.pb 

pressure and buffer pressure (Pa) 

- 

dimensionless quantity 

r 

offset distance from the crank to the center of shaft (m) 



R 

gas constant R = 8.314 (J/K mol) 

Subscripts 

R 

radius of gear (m) 

c 

compression chamber 

S 

stroke (m) 

d 

displacer 

T 

temperature (I<) 

e 

expansion chamber 

V 

volume (m 3 ) 

h 

heater 

W[, Wf, W s dimensionless indicated work, forced work, and shaft 

lc 

cooler 


work 

1 

lower yoke 

yo, yi 

lengths of displacer and piston rods (m) 

max, 

minmaximum and minimum 

yd.y P 

locations of top and bottom links of rhomboid (m) 

p 

piston 



r 

regenerator 

Greek symbols 

s 

swept volume 

«d 

displacement phase angle (deg) 

t 

total 

a v 

volume phase angle (deg) 

u 

upper yoke 

T 

compression ratio 




the )3 type possesses higher power density than the other two 
types. Therefore, the attention of this study is focused on the beta- 
type Stirling engine. Furthermore, a mechanism is typically 
installed with the beta-type Stirling engine to drive a displacer and 
a piston in order to fulfill the required phase angle between the two 
moving parts as well as the designed strokes of the displacer and 
the piston. The drive mechanism with the beta-type Stirling engine 
could be a mechanism of crank drive, bell-crank drive, twin inclined 
Scotch yoke drive, spring drive, or rhombic drive [ 6 , 24 — 26 ]. As was 


described in Ref. [ 26 ], the rhombic drive uses a jointed rhomboid to 
convert linear motion of the reciprocating piston and displacer to a 
rotation of the flywheel. In the rhombic drive mechanism, one rigid 
rod is installed connecting the piston to the top link of the jointed 
rhomboid, and another rigid rod connecting the displacer to the 
bottom link of the rhomboid. Two symmetric gears of equal 
diameter are connected to the right and the left links of the 
rhomboid fixed on the gears at an offset distance from respective 
gear centers. The two gears are in contact and rotate in opposite 


expansion 

chamber 




Fig. 1. Schematic and photo of 300-W beta-type Stirling engine with rhombic-drive mechanism [Model NCKU300-2, developed by Power Engines and Clean Energy Laboratory 
(PEACE Lab.), National Cheng Kung University]. 
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directions. The rhombic drive mechanism is a symmetric mecha¬ 
nism that makes the co-axial motion of the displacer and the piston 
move smoothly to minimize the lateral force acting perpendicular 
to the axis, and hence, the friction at contact surface between the 
piston and the cylinder can be greatly reduced. For this reason, the 
rhombic drive mechanism is frequently adopted by the (3-type 
Stirling engine designers. 

More recently, the present authors built a 300-W |3-type Stirling 
engine with a rhombic drive mechanism (Model NCKU300-2), of 
which the photograph is shown in Fig 1. According to the experiences 
learned from the development of the engine, it is recognized that 
major fundamental problems that the designers of the Stirling en¬ 
gines encountered are the time-consuming design process as well as 
the costly manufacturing procedure. In most of the cases the deter¬ 
mination of the geometric parameters of the rhombic drive is based 
on the experience of the designers or on insufficient experimental 
information. This is why the engine design can only be carried out in a 
relatively high-cost time-consuming trial-and-error process. Unfor¬ 
tunately, there are no desired explicit guidelines for design of the 
geometry of the rhombic drive mechanism in the existing reports. 

Under these circumstances, optimization of the rhombic drive 
mechanism used in the (S-type Stirling engine is attempted in this 
study based on a dimensionless theoretical model. Schmidt anal¬ 
ysis is incorporated with Senft’s shaft work theory [13,18] to build 
the dimensionless thermodynamic model, which is capable of 
providing the dimensionless shaft work under various combina¬ 
tions of the influential geometric parameters. The dimensionless 
thermodynamic model is verified with experimental data, and the 
optimal geometric parameters are so determined that the shaft 
work can be maximized. Eventually, the design charts that can be 
used as the guidelines for design of the geometry of the rhombic 
drive mechanism are provided. 

2. Dimensionless analytical model 

2.1. Displacements of piston and displacer 


e p - 1 > tp > 1 and e d - 1 > e d > 1 (3) 

Fig. 2(a) shows the snapshots of the typical rhombic drive 
mechanism in a cycle, which meets the above valid ranges of the 
parameters. However, if one of the criteria is not obeyed, a collision 
between the parts of the mechanism may occur. As firstly shown in 
Fig. 2(b), for the cases with e p < 1 (R — l u shorter than r) or f d < 1 
(R - l\ shorter than r), joint A is very possible to collide with the 
displacer rod at a high rotation speed. In Fig. 2(b), it is at 6 = 180° 
that the collision may take place. Secondly, as e p - 1 < e p (/ p - r 
smaller than R - / u ) or e d - 1 < e d (i d - r smaller than R - /i), the 
mechanism is not able to complete one cycle due to interference 
among the links (see 6 = 0°). 

Based on Eqs. (la) and (lb), the maximum and minimum 
heights of the top and the bottom links can be derived. The results 
are 


= \/(e„ + l) 2 - 


Tprnax 


yp.min ~ V ( e p 1) ' 


Td.max — V ( e d 1) e d 



(4a) 

(4b) 

(4c) 


(a) e-l>e >1 and 



In the modeling, it is assumed that (1) the transient variations in 
the volumes of the working spaces inside the engine are not si¬ 
nusoidal functions of time, but determined from the geometric 
relations among the parts of the mechanism, (2) the temperatures 
in the working spaces are unequal but constant, and (3) the pres¬ 
sures in the working spaces are varying with time periodically but 
are uniform at any time instant. 

The dimension variables and the schematic of the [3-type Stirling 
engine with the rhombic drive mechanism are shown in Fig. 1. 
According to the relations among the variables with the rhombic 
drive mechanism, the heights of top link of the rhomboid con¬ 
nected to the piston (y p ) and bottom link of the rhomboid con¬ 
nected the displacer (y d ) can be expressed in dimensionless form, 
respectively, as 

y p (0) = sin 6 + \je\ - (e p + cos d) 2 (la) 


y d (0) = sin 6 - \je\ - (e d + cos B) 2 


(lb) 


where the dimensionless parameters are defined as below: 

l n In R — hi R — /] 


v -^vj--e — — en — — e - R ~ lu 

y p y->y d y ' 


£ d 


(2a -f) 



e=w 


9 = 90" 


9 = 1 80” 


9 = 270" 






(c) e p -\<E p or e d -\<e d (invalid) 


It is noted that in practices the relations among these above 
dimensionless parameters must meet the following criteria 


Fig. 2. Valid and invalid parameter ranges, (a) e p - 1 > e p > 1 and - 1 > ed > 1 
(valid) (b) e p < 1 or e d < 1 (invalid) (c) e p - 1 < e p or e d - 1 < e d (invalid). 
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iM.min — \j ( e d + 1) ^ (4d) 

The dimensionless strokes with the piston and the displacer can 
be yielded as 

-> P = yp.max — yp.min (5) 

3 d = ^d,max — ^d,min (6) 

where S p = S p /r and S d = S d /r. Thus, one obtains the relation for 
the swept volume ratio 

K _ ^p = \Z( e p + 1) ~ f p ~ \/( e P ~ j) ~ e p ^ 7 j 

\J ( e d + 1 ) 2 “ e d “ V( e d - l ) 2 - f d 

The displacer leads the piston by a phase angle, which is readily 
determined by the difference in angle between y d max and y p max 
as 


a d = (180 /tc) 


Tt - cos 1 




- cos 1 



(in deg) 

( 8 ) 


Note that the height of the bottom surface of the displacer can 
be determined by y d + y 0 , and that of the top surface of the piston 
by y p + y 1 , where y 0 and Ja, are the dimensionless lengths of the 
displacer rod and the piston rod, respectively. The mechanism 
should be properly design to minimize the dead volume, and hence 
the minimum difference in a cycle between y d (0)+y o and 
y p (6) + y! is set to be zero such that the relation between y 0 andjq 
must satisfy the following equation 


yo -yi = \/( e P + e d ) 2 - (f P - e d ) 2 (9) 

Fig. 3 shows the dimensionless displacements of the bottom 
surface of displacer and the top surface of piston at e p = e d = 2.5 and 
e p = c d = 1. The crank angle at which the minimum difference in a 
cycle betweeny d (0) + yo andy p (0) + yi approaches zero is 180° or 
540° in this figure. It is also noted that the displacement phase 
angle between the displacer and the piston a d is approximately 
58.412°, which is obtained by using Eq. (8). 



Fig. 3. Dimensionless displacements of top surface of piston and bottom surface of 
displacer, at e p = e d = 2.5 and e p = e d = 1. 

Ve(6) = [yd,max-yd(0)1 Ad , 

L _ (11) 
Vciff) = {[yd(0)+yo]- [y P (0)+yi]}/s d 

where V e (8) = V e /V se and V c (8) = V c /V se . 

Then, dimensionless total volume of all working spaces in the 
engine can be written as 

F t (0) = v h +v I+ v k +{ (y d ,max+y 0 ) - [y p (0) +yi]}/s d ( 12 ) 

The first three terms represent the dimensionless dead volume 
in the heater, regenerator, and cooler, respectively. Note that the 
total volume variation is caused by the motion of the piston only. 
The ratio of the maximum to the minimum values of the total 
volume V t (8) denotes the compression ratio (y). Introducing a total 
dead volume ratiox = Vh + ^r + ^i; t° Eq. (12) yields 

y = 1 +KS d /[xS d + (y d ,max +y 0 ) - (y P .max +yi)] (13) 


2.3. Pressure and buffer pressure 


2.2. Volume variations of working spaces 

The volumes in the expansion space and the compression space 
caused by the movement of the piston and the displacer in an ideal 
situation with zero dead volume in the expansion space and the 
compression space can be written as 

14(0) = ^ d r[y d -y d (0)l , 

_ _ r_ _ i •) (10) 

Vc(0) = ^ P r{[y d (0) +y 0 ] - [y P (0) +yij} 

where A d and A p are the cross section areas of the displacer and the 
piston, respectively. To simplify the model, since the cross-sectional 
area of the rod connected to the displacer is much less than that of 
the displacer (about 0.5 ~ 2%), one may ignore the cross-sectional 
area of the rod and let A p ==A d . Also, let \4e be the swept volume by 
the displacer (= A d rS d ). One can derive the dimensionless volume 
variations in the expansion space and the compression space with 
zero dead volume as 


In this study, Schmidt theory is employed to determine the 
pressure in the working spaces as 


P(0) 


mR 


Vh , v r V k V'(6) V c (8) 

T^TT-TT-I T-1 T- 

ih if lk l e 1c 


(14) 


where T e and T c denote the temperatures in the expansion space 
and the compression space, respectively. In addition, T k , T r , and T k 
are the temperatures in the inevitable dead volumes inside the 
heater, the regenerator, and the cooler. The temperatures in the 
dead volumes may be further approximated by the mean value of T e 
and T c . Eq. (14) can then be simplified and written in dimensionless 
form as 


p(6) = p(0) = -C—-=- (15) 

/ Vse 2l X /0+T) + TV e (d) + V c (d) 

where t = T c /T e is referred to as the temperature ratio and p(8) is the 
dimensionless pressure. As described by Cheng and Yang [27], the 
temperature ratio i is a function of thermal characteristics of the 
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heater and the cooler, effectiveness of the regenerator, mass of the 
working gas, and engine speed. 

The dimensionless buffer pressure p b in the crank case (or buffer 
space) needs to be determined. Assume the volume in the crank 
case is sufficiently large; therefore, the pressure variation caused by 
the piston motion can be neglected. Thus, the dimensionless buffer 
pressure p b can be approximated as a constant. Herein, the 
dimensionless buffer pressure is given by the geometric mean of 
Pmax ar, d Pinin' \/PmaxPmin- 


2.4. Indicated, forced, and shaft works 


Dimensionless indicated work, forced work and shaft work are 
defined as: 


Wi = 


Wj 

mRT e 


W f 


W f 

mRT e 


W s 


W s 

mRTe 


(16) 


With the help of Eqs. (12) and (15), the p - V t diagram for the 
case with x = 0.3, % = 0.5, e p = e b = 2.5 and e p = e d = 1 are obtained 
and plotted in Fig. 4. The dimensionless buffer pressure p b is indi¬ 
cated in the figure. The area inside the p - V t is evaluated by 
integration fpdV t to yield the dimensionless indicated work W- v 




(17) 


Forced work is the work loss due to work transmission via the 
mechanism. In a fl-type Stirling engine, forced work only acts on the 
piston and is only dependent on the motion of the piston, resulting 
from the pressure difference on the both sides of the piston. Based 
on the definition given by Senft [ 28 ], the dimensionless forced 
work can be calculated by. 


W f 


P - Pb d T] 


dd 




(18) 


where the symbol [] means 


{ -a as a < 0 
| 0 as a > 0 


(19) 


Finally, as described in Ref. [ 28 ], one can determine the 
dimensionless shaft work by 


W s = EWj -(1/E- E)W { (20) 

where £ is the effectiveness of the mechanism, depending on in¬ 
dividual mechanism design, friction, engine speed, lubrication, and 
so on. The magnitude of £ is usually determined by experiments 
whose typical value is ranged from 0.7 to 0.9. In this study, the value 
of £ is assigned to be 0.8. 

3. Experiment and comparison 

In this study, performance measurement of Model NCKU300-2 
is conducted, and the experimental results are utilized to partially 
compare with theoretical results for validation. The basic specifi¬ 
cations of Model NCKU300-2 are tabulated in Table 1, and the 
design parameters of the links are shown in Fig. 5 . More detailed 
information of the experiment apparatus has been described in 
Ref. [ 29 ], The engine is heated by a 1.21-kW infrared heater, and a K- 
type thermocouple and a PID controller are equipped with the 
heater in measuring and monitoring the heating temperature, 
respectively. A torque sensor and an angular speed transducer are 
installed in between the engine and a hysteresis brake. As the data 
of the torque and the angular speed are recorded, shaft power of 



Fig. 4. Dimensionless pressure variations and buffer pressure versus dimensionless 
total volume variations, at t = 0.3, % = 0.5, e p = e d = 2.5 and e p = e d = 1. 


Table 1 

Specifications of model NCKU300-2. 


Parameter 

Value 

Power capacity 

300-500 W 

Stroke, S d 

35.219 mm 

Cross section area of displacer, Ad 

38.485 cm 2 

Displacer swept volume, V se 

135 c.c. 

Displacement phase angle, a d 

54.219 

Volume phase angle, a v 

105.537 

Swept volume ratio, k 

1 

Dead volume ratio, x 

1.478 

Compression ratio, y 

1.364 

Effectiveness of mechanism, E 

0.73 

Buffer pressure, p b 

8 atm 

Mass of working fluid 

0.0752 mol 

Working fluid 

Helium 

Heater type 

Tubular heater 

Heating temperature 

900-1100 K 

Cooler type 

Fined water jacket 

Cooling temperature 

300 K (liquid cooling) 
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Fig. 5. Comparison between theoretical and experimental results. 

the engine can be determine from the product of the torque and the 
angular speed. As described in Refs. [ 27 , 29 ], at specific heating and 
cooling temperatures, there exists a maximum shaft power at a 
certain engine speed. The maximum shaft power and the temper¬ 
ature ratio under different heating temperatures are measured by 
experiments. Finally, the plot of dimensionless shaft work W s 
versus temperature ratio z is yielded and shown in Fig. 5 . In this 
figure, it is found that the temperature ratio is located in the range 
from 0.49 to 0.58. It is noted that the experimental data agrees with 
the theoretical results with a relative error lower than 5.2%. The 
feasibility of the theoretical model is verified; therefore, in the 
following section a parameter analysis is performed and the geo¬ 
metric optimization of the engine is attempted. 

It is noted that the effectiveness of mechanism E can be varied 
between 0.7 and 0.9 typically. In the pure theoretical analysis sec¬ 
tion of this study, the value of E is chosen to be the mean value, 0.8 
simply for demonstrating a standard parametric analysis. However, 
as the theoretical model is applied for a specific engine, this value 
should be determined by fitting the experimental data exclusively 



Fig. 6. Effects of temperature ratio on dimensionless works, at E = 0.8, x = 0-5, 
= = 2.5, and £ p = e d = 1. 



Fig. 7. Effects of total dead volume ratio on dimensionless works, at E = 0.8, x = 0.3, 
e p = e d = 2.5, and e p = e d = 1. 



(a) Effects of e p 



Fig. 8. Effects of e p and ed on dimensionless works, at E = 0.8, x = 0.3, % = 0.5, and 
e p = e d = 2.5. 
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W 


0.14 
0.12 
0.1 
0.0 8 


(a) £ = 0.7 




(c) E = 0.9 

Fig. 9. Dimensionless shaft work as a function of e p and e A for different effectiveness of 
mechanism at t = 0.3, % = 0.5, and c p = td = 1. (a) E = 0.7, (b) £ = 0.8, (c) E = 0.9. 

from the performance experiments conducted for the engine. 
Table 1 shows the specifications of Model NCKU300-2, whose 
effectiveness of mechanism £ is found to be 0.73 in accordance with 
the experiments. 

4. Results and discussion 

According to the discussion in the precedent sections, the 
dimensionless shaft work W s is a function of E, r, x, e p , e d , e p , and c d . 
As has been described earlier, the magnitude of E is determined by 


(a) 7 = 0.5, J = 0.5 


0.12 
w . 0.1 
0.08 


(b) 7 = 0.3, x = 1° 

Fig. 10. Dimensionless shaft work as a function of e p and e d for different combinations 
of r and x. at £ = 0.8 and e p = e d = 1. (a) x = 0.5, % = 0.5, (b) x = 0.3, % = 1-0. 

experiments, whose typical value is ranged from 0.7 to 0.9. Here in 
this study, the value of E is determined to be 0.8. The present study 
is aimed to seek optimal combination of these influential param¬ 
eters toward maximization of the dimensionless shaft work W s . 

Fig. 6 shows the effects of the temperature ratio x on W- v W s , and 
Wf, at E = 0.8, x = 0.5, e p = e d = 2.5, and e p = £ d = 1. It is found that 
both W j and W s are decreased, while W f is increased, by increasing 
the temperature ratio. As r is elevated to be 0.81, the magnitude of 
W s vanishes to zero. This implies that the engine is not able to 
output and work. As the temperature ratio is further increased to be 
unity, even the indicated work will become zero since x = 1 means 
no temperature difference between the heater and the cooler. 

Effects of the total dead volume ratio x on VV;, W s , and W f are 
displayed in Fig. 7 for the case at £ = 0.8, x = 0.3, e p = e d = 2.5, and 
£p = £d = 1. The results shown in this figure clearly reveal that the 
works produced by the engine are sensitive to an increase in the 
dead volume. It is observed that W,, W s , and W f are all remarkably 
decreased by increasing the dead volume ratio x . especially in the 
regime between x = 0 and 0.5. In theory, the magnitude of x should 
be reduced as low as possible. However, installation of the heater, 
the regenerator, and the cooler is necessary in an engine in order to 
maintain a temperature difference (and hence, decrease the 
magnitude of t). These three components and their internal dead 
volumes are actually inevitable, and in this regard, a study of the 
optimal value of the total dead volume ratio x becomes an impor¬ 
tant issue. 

Effects of e p and £ d on the dimensionless works are conveyed in 
Fig. 8(a) and (b) respectively, for the case at £ = 0.8, i = 0.3, x = 0-5, 
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T 


(a) Optimal e d 



T 


(b) Optimal e p 



T 


(c) Maximum IV 

Fig. 11. Dependence of optimal W s , e p , and e d at E = 0.8 and e p = e d = 1. (a) Optimal e d , 
(b) Optimal e p , (c) Maximum.W s . 


and e p = e d = 2.5. In Fig. 8(a), it is found that W, and W f are not 
sensitive to a change in e p . As e p is elevated from 1.0 to 1.25, either 
Wj or W { experiences an inappreciable increase. However, a slight 
decrease in W s is observed as e p is elevated from 1.0 to 1.25. On the 
contrary, Fig. 8(b) shows that both W, and W { are decreased by 
increasing and Wj is decreased faster than W f . As a result, W s is 
decreased as ed is increased. The results in Fig. 8(a) and (b) show 
that all the works are monotonically varied with either e p or ed. It 
appears that both e p and ed have to be set as small as possible to 
yield a high shaft work. Being aware that in accordance with the 
criteria of Eq. (3), e p and ed ought to be higher than or equal to unity. 
Therefore, the optimal values of e p and ed both are suggested to be 
1 . 0 . 

Fig. 9 displays three-dimensional surface of W s as a function of 
e p and ed for different £, at i = 0.3, % = 0.5, and e p = ed = 1. It is 
clearly seen in Fig. 9(a) and (b), for the cases at £ = 0.7 and 0.8, that 
a peak of W s exists in each of the plots. In Fig. 9(a), the optimal 
combination of e p and ed is e p = 2.24 and ed = 3.32. However, in 
Fig. 9(c) for the case with a higher £, the peak of W s is no longer 
visible. In this case, W s is increased monotonically by a decrease in 
e p . Therefore, under the constraint of Eq. (3), the optimal e p is 
selected to be 2.0. On the other hand, at e p = 2, it is observed in 
Fig. 9(c) that W s increases monotonically with e p . The selection of 
ed is dependent on individual cases. However, a finite value ought 
to be assigned to e p due to the limitation of links size. In this study, 
ed is suggested to be 5.0. 

Fig. 10 conveys the results of W s as a function of e p and ed for 
different combinations of x and x, at £ = 0.8 and t p = f d = l.A closer 
look into Fig. 9(b) and Fig. 10(a) gives that as the value of t is 
assigned to be 0.3 or 0.5, a peak of W s is clearly seen. The peak is 
more obvious with a higher t. Fig. 10(b) exhibits a peak of W s at a 
small e p . A comparison between Fig. 9(b) and Fig. 10(b) illustrates 
the influence of x varied from 0.5 to 1.0. The points of the peaks of 
W s indicate the optimal combinations of e p and ed under different 
conditions. 

Next, the optimal values of e p , ed, and the corresponding 
maximum W s are plotted as design charts that may be used to 
determine optimal geometric parameters and evaluate the shaft 
power output. 

Dependence of the optimal values of e p , ed and their corre¬ 
sponding W s on t andx are displayed in Fig. 11 at £ = 0.8 and 
e d = e p = 1.0. Using the plots of this figure, one is able to find out the 
optimal values of e p , e d and W s provided that the values of i, and x 
are given beforehand. For example, if i = 0.5 and x = 0.5, the 
optimal design of rhombic drive mechanism is: e p = 2.12 and 
e d = 2.85, and the corresponding W s is 0.15. It is noted that for the 
curve ofx = 0.1, e d is fixed at 5.0 when t is smaller than 0.42. 

5. Conclusions 

A dimensionless model has been developed and used to opti¬ 
mize the geometry of the rhombic drive mechanism with the (3- 
type Stirling engine toward maximization of the shaft work. Dis¬ 
placements of the piston and the displacer with the rhombic drive 
mechanism and variations of volumes and pressure in the cham¬ 
bers of the engine are expressed in dimensionless form. Schmidt 
analysis is incorporated with Senft’s shaft work theory to build the 
dimensionless thermodynamic model, in which the dimensionless 
shaft work W s is treated as a function of the effectiveness of 
mechanism (£), the temperature ratio (t), the total dead volume 
ratio (x), and the geometric parameters (e p , e d , e p , and e d ). In order 
to verify the thermodynamic model, the performance test of a 300- 
W (3-type Stirling engine with rhombic drive mechanism is also 
conducted. It is noted that the experimental data agrees with the 
theoretical results with a relative error lower than 5.2%. 
















978 


C.-H. Cheng, H.-S. Yang /Energy 64 (2014) 970-978 


In addition, a parameter analysis is performed and the geo¬ 
metric optimization of the engine is attempted. The design charts 
that help determine the optimal geometric parameters of the 
rhombic drive mechanism are provided. 

Results for a typical case with E = 0.8, % = 0.5, e p = e d = 2.5, 
and e p = f d = 1 show that as the temperature ratio i is elevated to 
be 0.81, the magnitude of dimensionless shaft work output W s 
vanishes to zero. In addition, the works produced by the engine are 
sensitive to an increase in the dead volume. It is observed that W- x , 
W s , and W f are all remarkably decreased by increasing the total 
dead volume ratio x. especially in the regime between x = 0 and 
0.5. 

Valid ranges of the geometric parameters of the rhombic-drive 
mechanism are found to be e p - 1 > e p > 1 and e d — 1 > e d > 1. 
In accordance with the predictions by the model, the optimal values 
of e p and £ d both are selected to be 1.0. On the other hand, the 
optimal values of e p , e d , and the corresponding maximum W s are 
plotted as design charts that may be used to determine optimal 
geometric parameters and evaluate the shaft power output. For the 
particular case at t = 0.5 and x = 0.5, the optimal combination of e p 
and e d is: e p = 2.12 and e d = 2.85, and the corresponding W s is 0.15. 
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